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Two-dimensional (2D) equations describing the nonlinear interaction between upper-hybrid and 
dispersive magnetosonic waves are presented. Nonlocal nonlinearity in the equations results in the 
possibility of existence of stable 2D nonlinear structures. A rigorous proof of the absence of collapse 
in the model is given. We have found numerically different types of nonlinear localized structures 
such as fundamental solitons, radially symmetric vortices, nonrotating multisolitons (two-hump 
solitons, dipoles and quadrupoles) , and rotating multisolitons (azimuthons) . By direct numerical 
simulations we show that 2D fundamental solitons with negative hamiltonian are stable. 
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I. INTRODUCTION 

Upper-hybrid (UH) waves are frequently observed in 
space and laboratory plasmas. The UH waves can be ex- 
cited by beam instabilities, mode conversion of extraor- 
dinary electromagnetic waves at the upper-hybrid reso- 
nance layer, etc J* One-dimensional (ID) theory of the 
nonlinear UH waves interacting with the low-frequency 
motions of magnetosonic type was developed in Refs. 
2-7. In particular, for low frequency magnetohydro- 
dynamic perturbations with frozen-in field lines and 
for the negative group dispersive UH waves, Kauf- 
man and Stenflo 2 - showed the existence of UH soli- 
tons with compressional density (magnetic field) pertur- 
bations in the super-magnetosonic regime. Dispersive 
magnetosonic/lower- hybrid waves interacting with high- 
frequency one-dimensional UH waves were first consid- 
ered in Refs. 4 and 5 and then were studied in more 
detail, including intensive numerical modelling 7 . 

The aim of this paper is to present two-dimensional 
(2D) theory of nonlinearly coupled dispersive magne- 
tosonic and high-frequency UH waves. We have derived 
a set of corresponding nonlinear equations and found 
numerically different types of 2D nonlinear localized 
structures such as fundamental solitons, radially sym- 
metric vortices, nonrotating multisolitons (dipoles and 
quadrupoles), and rotating multisolitons (azimuthons). 

Dispersion of the magnetosonic wave effectively intro- 
duces a nonlocal nonlinear interaction, that is the non- 
linear response depends on the wave packet intensity at 
some extensive spatial domain. Nonlocal nonlinearity 
naturally arises in many areas of nonlinear physics and 
plays a crucial role in the dynamics of nonlinear coher- 
ent structures. While collapse is a usual phenomenon 
in the multidimensional Zakharov-like models with lo- 
cal low-frequency responce, nonlocal nonlinearity can ar- 
rest collapse and lead to stable multidimensional local- 
ized structures ^ 10 ' 11 : 12 ' 13 . Turitsyn proved-^ the ab- 
sence of collapse for three particular shapes of the non- 
local nonlinear response in the multidimensional gener- 
alized nonlinear Schrodinger equation (GNSE). Later, a 
rigorous proof of absence of collapse in arbitrary spa- 



tial dimensions during the wave-packet propagation de- 
scribed by the nonlocal GNSE with sufficiently general 
symmetric response kernel was presented in Ref. 15. In 
the present paper we give a proof of absence of collapse 
for the 2D model describing the nonlinear interaction be- 
tween upper- hybrid and dispersive magnetosonic waves. 

The paper is organized as follows. In Sec. [Til we present 
the generalized Zakharov-type system of 2D equations, 
describing the interaction between high-frequency upper- 
hybrid waves and low-frequency dispersive magnetosonic 
waves. A linear stability analysis is performed in Sec. IIIII 
Section IIVI contains the rigorous proof of absence of col- 
lapse in the model. Localized nonlinear solutions, includ- 
ing vortex ring-like, nonrotating (a monopole, a dipolc, 
and a quadrupole) and rotating (azimuthons) multisoli- 
ton solutions are presented in Sec. [V] The conclusion is 
made in Sec. I VII 



II. DERIVATION OF EQUATIONS 

We consider a homogeneous electron-ion plasma in a 
uniform external magnetic field Bo = Bqz, where z is the 
unit vector along the z direction. In the linear approxi- 
mation, the UH waves are characterized by the dispersion 
relation 
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Equation for the slow varying complex amplitude ip of 
the potential of the high-frequency electrostatic electric 
field 
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of the upper hybrid wave can be obtained from the equa- 
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where the dielectric tensor e is considered as a differential 
operator with lu — * u>uh + id/dt (assuming louh d/dt) 
and k — * — iV. Under this, the nonlinear perturbations 
of the plasma density Sn and magnetic field SB are taken 
into account in e±, eii and g, so that substitutions no — > 
ri + 5n and £? — > -Bo + <5.B are made in the resulting 
equation 
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The second term in the {. . . } bracket in Eq. © comes 
from the last term in Eq. © and corresponds to the so- 
called vector nonlinearity. This term is identically zero 
for ID case and for the fields with axial symmetry. This 
term can also be neglected if uj ce /ujp e <C 1. 

The upper-hybrid waves have wave numbers almost 
normal to the external magnetic field (k z <C k±) and, 
in the following, we will consider two-dimensional (2D) 
case with k z = so that A = d 2 /dx 2 + d 2 /dy 2 and 
V = (d/dx,d/dy). 

The low-frequency motion of the plasma is governed 
by the continuity and momentum equations for ions and 
electrons. We assume quasineutrality condition, so that 
Sri; = 8n e = Sn. Thus, we have 
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where ji (j e ) is the ion (electron) ratio of specific heats, 
and 
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are nonlinear terms in electron equations, the angular 
brackets denote averaging over the high-frequency os- 
cillations, and the superscripts H denote corresponding 
quantities for the high-frequency fields. Multiplying Eq. 
(flTjf by m/M and adding with Eq. (jHJ), we have 
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where we have introduced the effective sound speed v s = 
y/{jiTi + j e T e )/M. Taking the div from Eq. JEJ) we get 
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Using Eqs. (J7J) and ©, and eliminating V x (v e — v$) 
with the aid of the Maxwell equation (we neglect the 
displacement current for the low-frequency motion) 
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where va = Bo/^/iimoM is the Alfven speed. The sec- 
ond term in the right-hand side of Eq. (| 1 5[) is small com- 
pared to the first one by the factor ~ lo/lojjjj and can 
be neglected. To obtain equation for the low-frequency 
magnetic field perturbation SB we substract Eq. (|10|) 
from Eq. ((SJ) and take the curl 
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Using Eqs. 0,©, (HI and the Maxwell equation 
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from the high-frequency momentum equation for elec- 
trons we have 
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With the aid of the Maxwell equation for 9B ff /3t, two 
terms in the expression for Fi can be combined to yield 

Fi = ((vf -V).vf + [vf x[Vxvf]]) 

= ^(V(vf -vf)). (21) 

Using Eqs. (| 1 5[> . (fT8"|) - (f2"Tj) we obtain equations for the 
low-frequency plasma density and magnetic field pertur- 
bations 
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In the linear approximation, Eqs. (f2"2"| and (f2"3")l give the 
dispersion relation for the dispersive fast magnetosonic 
wave 
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In what follows we omit the subscript _L in k± . Equations 
© , (|2"2"|) and (|2"3")) form a closed system of 2D equations 
describing the interaction between upper-hybrid waves 
and dispersive magnetosonic waves. In the ID case these 
equations coincide with those obtained in Refs. 4, 5 and 
7. 



III. NONLINEAR DISPERSION RELATION 

In this section we consider the linear theory of the mod- 
ulational instability of a pump wave with a frequency 
close to the upper-hybrid frequency. As usual, we express 
the low-frequency perturbation of the plasma density as 



= ncxp(ik • r — iilt) + c.c. 
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while the upper-hybrid wave is decomposed into the 
pump wave and two sidebands 

+^ e «K k »- k )- r -(*»- !! )'l+c.c, (27) 



where Sq = louh k 2 R 2 /2. The amplitudes of the satellites 
can be calculated from Eq. ([B]). We have 
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and the function 



D ± = 2u; UH (k ±k) 2 {S ± T ^) 



(30) 
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is the Fourier transform of the linear operator in the left- 
hand side of Eq. © evaluated at k ± k, and S± = 
W[/_ffi? 2 [(ko ± k) 2 — /cq]/2 are the mismatches between 
the satellite frequencies and the frequency of the pump 
wave. The amplitudes of the low-frequency perturbations 
are found from Eqs. ((22j) and ([23jl : 
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and il 2 is determined by Eq. (f24|) . By combining Eqs. 
([28|) . (f29|) and (|32|) one can obtain a nonlinear dispersion 
relation 
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Equation (|35|) generalizes the nonlinear dispersion rela- 
tion obtained by Eliasson and Shukla 7 by considering 2D 
case and including the vector nonlinearity. For ID case 
we recover the previous result. Note, that in the case 
of coplanar (in the plane perpendicular to the magnetic 
field) wave vectors k || ko, the parametric coupling of 
the waves due to the vector nonlinearity is absent, while 
the coupling due to the scalar nonlinearity is the most 
effective. In the opposite case, i.e. k _L ko, the interac- 
tion due to the vector nonlinearity is the most effective, 
while the interaction due to the scalar nonlinearity is the 
least effective (and absent for fco <C fc). In general case 
k D" ko, and when oj pe ~ uj ce so that both types of the 
nonlinearities yield comparable contribution, this leads 
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to a rather complicated picture of the parametric insta- 
bility. In this paper, we restrict ourselves to the case of 
a weakly magnetized plasma with uj 2 e /u 2 e <C 1. Results, 
concerning the case of a moderately magnetized plasma 



with i 



u pe, 



including the nonlinear analysis (see be- 



low), will be published elsewhere. 

In the case uj ce <C uj pe , one can neglect the vector non- 
linearity and Eq. (|35p takes the form 
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When both sidebands are resonant and the low-frequency 
perturbations are nonresonant, we can consider the lim- 
iting case ko 3> k and Q/. 3> f2. Then, the dispersion 
relation Eq. I|36p takes the form 
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where v g is the group velocity of the UH wave. Equation 



([38]) predicts the instability when |^ | > 8Trn MR 2 nf.. 
In the opposite case of long-wavelength pump k ^ ko 
Eq. ([36]) is reduced to 
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where 5 = ujjjh k 2 R 2 /2 and, thus, there is a purely grow- 
ing instability for sufficiently large amplitudes of the 
pump wave Eq. 



IV. PROOF OF ABSENCE OF COLLAPSE 

In what follows, we will consider the case ui 2 e /uj 2 e <C 1, 
so that vector nonlinearities in Eqs. © and (|22[) can be 
neglected. Then, introducing the dimensionless variables 
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and eliminating Sn/n a , one can rewrite Eqs. ([5]), (122[) 
and (|23p as follows 
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Equations (|43|) and (|44f conserves energy 

N = I |V^| 2 dr, (45) 



and Hamiltonian 
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In this section, following the ideas suggested in Ref. 
14, we present a rigorous proof of the absence of collapse 
for the stationary 2D solutions of the form tp(x,y,t) = 
ip(x, y) exp(iAi) in the model described by Eqs. (|43|) and 
(|H]) . We use an exact approach 1 ^ based on the Liapunov 
stability theory. Let us briefly recall the essence of the 
Liapunov method. For an invariant set of the dynami- 
cal system (in particular, for the set of stationary solu- 
tions {u s }) to be stable, it is sufficient that there exists a 
functional L[u] with the following properties: a) L is the 
positive definite functional for the perturbed states {«}, 
i.e. L[u] ^ 0; b) L reaches its minimum on the set {u s }, 
L[u s ] = 0; c) L is nonincreasing function of time t, i.e. 
dL/dt ^ 0. For hamiltonian systems, these conditions 
arc equivalent to the requirement that the Hamiltonian 
is bounded from below under the fixed conserved quan- 
tity N (then, one can choose L = H — i? m , n ). Then, the 
2D stationary localized solutions corresponding to the 
global minimum of H (fundamental solitons) are stable 
in the Liapunov sense. 

Note, that for the stationary solutions, Eqs. (|43[) and 
can rewritten as follows 



A(-A/i + A)tp = V • [(/3 - A)bVip], 



with 
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where G(r) — Ko(y/(3 + l|r|)/27r is the Green function 
satisfying the equation (/3+1 — A)G(r) = <5(r), and K (z) 
is the modified Bessel function of the second kind of or- 
der zero. It is seen, that the nonlinearity in Eq. (|48|) has 
essentially nonlocal character. We represent the Hamil- 
tonian (l46l) as 



H = A 



1 



(50) 



where A = J \Aip\ 2 dr and the nonlinear term P can be 
written as 
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where 

C = J |V^(r)| 2 |V^(r')| 2 G(|r-r'|)drrfr' > 0. (52) 
To get an estimate for C we use the following inequality 14 

/^*< 2 (//**) 1/2 (/(v/)**) 1/2 , 

(53) 

where /(r) is an arbitrary sufficiently smooth function, 
the integration in Eq. (|53[) is performed over entire 2D 
space. We have 

C= I |V^(r)| 2 ^^G(|r-r1)|r-rVrdr' 
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^ M J |Vv?(r)| 2 dr J 
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= 2A1N 3/2 A 1/2 , (54) 
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where M = max 2 [zA'o(^\//? + 1 )/2i"] > and we have 
used Eq. ([53)) . Next we use the inequality 



|Vv3(r)| 4 dr < iV 2 . 



(55) 



As a result, one can write the obvious chain of inequali- 
ties: 

P > -N 2 + C ^ -N 2 -C ^ -N 2 -2MN 3 / 2 A 1 ' 2 . (56) 

Inserting this estimate into the expression for the hamil- 
tonian Eq. ([50]) we get 



(57) 



H ^ A - -N 2 - MN 3/2 A 1/2 . 
2 



Under the fixed N, the right hand side of the inequality 
([57]) reaches its minimum at A = M 2 N 3 /4, so that 



H ^ —M N /4 — N /2. 



(58) 



Thus, we have showed that, under the fixed conserved 
quantity N, the Hamiltonian is bounded from below. Ac- 
cording to standard Liapunov theory, this represents a 
rigorous proof that a collapse with the wave amplitude 
locally going to infinity cannot occur for the stationary 
solutions in the model described by Eqs. (|43|) and (1331). 



V. NONLINEAR LOCALIZED SOLUTIONS 

We consider the case when the upper-hybrid wave- 
length is much smaller than the typical spatial scale of 
the localization and represent (in the corresponding di- 
mensionless variables) 



We assume that koL ^> 1, where L = is the typ- 

ical length scale of the low-frequency perturbation. We 
also make an additional assumption of ko • V = 0. Under 
this, the group velocity of the upper-hybrid wave is per- 
pendicular to the direction of inhomogeneity and there is 
no transportation of energy in the direction of inhomo- 
geneity which results in stationary envelope structure in 
the direction perpendicular to the magnetic field. Then, 
the system Eqs. (|43|) and (|44f takes the form 



(60) 



+ A$ = (p - A)M>, 

ot 
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where we have rescaled $ so that $ — > ko$. The proof of 
the absence of collapse for the stationary solutions of the 
system (jSU)) and (pjlj) is quite analogues to that presented 
in the previous section. 

We look for stationary solutions of Eqs. (16T))) and (|6"Tj) 
in the form <$>(x,y,t) — ty(x,y) exp(i\t/[i), where A//z is 
the nonlinear frequency shift, so that obeys the equa- 
tion 



A* + A* = (/?- A)6tf, 



(/?+!- A)6= -|*| 2 



(62) 



(63) 



To solve numerically Eqs. (f6"2")) and (f6"3")) . we impose peri- 
odic boundary conditions on Cartesian grid and use the 
relaxation technique similar to one described in Ref. 17. 
Choosing an appropriate initial guess, one can find nu- 
merically with high accuracy (the norms of the residuals 
were less than 10 -9 ) three different classes of spatially 
localized solutions of Eqs. lf6"2"]) and (f6"3"|) - the nonrotat- 
ing (multi)solitons, the radially symmetric vortices, and 
the rotating multisolitons (azimuthons). 

The real (or containing only a constant complex fac- 
tor) function ^(x,y) corresponds to nonrotating solitary 
structures. Examples of such nonrotating (multi)solitons 
for Eqs. l|6"2")) and l|6"3"j). namely, a monopole, a dipole, 
two-hump soliton, and a quadrupole are presented in 
Figs.QTa)-[]Td) for the case /3 = 0.1. The analogous solu- 
tions can be found for (5 > 1. The nonrotating multipolcs 
consist of several fundamental solitons (monopoles) with 
opposite phases. 

The second class of solutions, vortex solutions, are 
the solutions with the radially symmetric amplitude 
|^(x, y)|, that vanishes at the center, and a rotating spi- 
ral phase in the form of a linear function of the polar an- 
gle 9, i.e. arg $ = m6, where m is an integer. The index 
m (topological charge) stands for a phase twist around 
the intensity ring. The important integral of motion asso- 
ciated with this type of solitary wave is the z-component 
of the angular momentum 



ip = $exp(ik ■ r — i<5 £). 



(59) 



M z = 



Im J [$*(r x V ± $)] z dr, 



(64) 
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FIG. 3: Examples of the rotating multisolitons (azimuthons) 
with A = 5 for /3 = 0.1: a) azimuthon with two intensity 
peaks; b) azimuthon with four intensity peaks. The field in- 
tensity 1^1 is shown. 



FIG. 1: Numerically found nonrotating stationary localized 
solutions of Eqs. (|60[1 and (|6ip for /3 = 0.1: (a) monopole 
with A = 0.5; (b) dipole with A = 0.5; (c) two-hump soliton 
with A = 0.5; (d) quadrupole with A = 2. The real part of 
the field ^ is shown. 
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FIG. 2: Examples of the stationary radially symmetric vortex 
solutions, characterized by the topological charges m = 1,2 
for A = 0.05 (left panel) and A = 1 (right panel): radial 
profiles of the field intensity (solid curve) and density n 
(dashed curve) are shown for f3 — 0.25. 



which can be expressed through the soliton amplitude U 
and phase 0, 



My 



d( t>TT2 



00 



U z dr, 



(65) 



and for the vortex we have M z — mN. Examples of the 
stationary radially symmetric vortex solutions, charac- 
terized by the topological charges m = 1, 2 for different 
values of A are shown in Fig. [2j 

The third class of solutions, rotating multisolitons with 
the spatially modulated phase, were first introduced in 
Ref. 18 for models with local nonlinearity, where they 



were called azimuthons. The azimuthons can be viewed 
as an intermediate kind of solutions between the rotating 
radially symmetric vortices and nonrotating multisoli- 
tons. Using variational analysis to describe azimuthons, 
the authors of Ref. 12 considered the following trial func- 
tion in polar coordinates (r,6) 

V(r,6) = r |m| $(r) (cos m9 + ip sinm0), (66) 

where $ is the real function, which vanish fast enough 
at infinity, m is an integer, and < p < 1. The case 
p = corresponds to the nonrotating multisolitons (e. 
g. m = 1 to a dipole, m — 2 to a quadrupole etc.), 
while the opposite case p = 1 corresponds to the radially 
symmetric vortices. The intermediate case < p < 1 
corresponds to the azimuthons. The trial function Eq. 
(|66p was chosen as an initial guess in our numerical re- 
laxation method. Then, the numerically found complex 
function ^(x,y) with a spatially modulated phase corre- 
sponds to the azimuthons. We introduced the parameter 
p (modulational depth), which is similar to the one in 
Eq. (f66|) . in the following way 



p = max |Im ^|/ max |Re^| 



(67) 



For fixed A, there is a family of azimuthons with dif- 
ferent p. Since the azimuthons have a nontrivial phase, 
they, like the radially symmetric vortices, carry out the 
nonzero angular momentum. In Figures [3] we demon- 
strate two numerically found examples of the azimuthons 
with two and four intensity peaks for the nonlocal model 
described by Eqs. (|6^|) and (f6"5| . The azimuthon with 
two intensity peaks consists of two dipole-shaped struc- 
tures in the real and imaginary parts of VP with different 
amplitudes and the ratio of these amplitudes is the mod- 
ulational depth p. Note, that the choice of initial guess 
in the relaxation method for finding the azimuthons is 
much more sophisticated than that for the nonrotating 
multisolitons or vortices. For example, we were not able 
to find azimuthon solutions with predetermined (in ad- 
vance) value of p. 

We next addressed the stability of these localized solu- 
tions and study the evolution of the solitons in the pres- 
ence of small initial perturbations. We have undertaken 
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extensive numerical modeling of Eqs. ((60)) and (|6"Tj) ini- 
tialized with our computed solutions with added gaus- 
sian noise. The initial condition was taken in the form 
^{x, y)[l+ef(x, y)}, where ^(x, y) is the numerically cal- 
culated exact solution, f(x, y) is the white gaussian noise 
with variance a 1 = 1 and the parameter of perturbation 
e = 0.005 -i- 0.01. In addition, azimuthal perturbation 
of the form is sin 9 was taken for the vortices and az- 
imuthons. Spatial discretization was based on the pseu- 
dospectral method. Temporal f-discretization included 
the split-step scheme. 

Numerical simulations clearly show that the funda- 
mental solitons are stable and do not collapse even for 
the negative initial hamiltonian. Stable evolution of the 
monopole soliton with A = 1 is shown in Fig. H{a). We 
have observed neither stable evolution nor collapse for 
multisolitons. If the nonlinear frequency shift A is not 
too large, the multisolitons decay into several monopole 
solitons, but can survive over quite considerable times. 
Splitting of the dipole soliton in two monopoles which 
move in the opposite directions without changing their 
shape is shown in Fig. |U(b). Figure 5 presents an ex- 
ample of the decay of the vortex into three fundamental 
solitons. Since the total angular momentum is conserved, 
the monopole solitons fly off the ring along tangential 
trajectories. A similar behavior was observed for the ro- 
tating multisolitons (azimuthons with two and four in- 
tensity peaks). If the nonlinear frequency shift A exceeds 
some critical value depending on the multisoliton type 
(i.e. dipole, azimuthon with two intensity peaks etc.), the 
unstable multisoliton turn into the one monopole with 
larger amplitude. Although the multisolitons turn out to 
be unstable in all our runs, one cannot exclude the possi- 
bility of existence of stable multisolitons in some narrow 
range of parameters (see, for example, Refs. 19 and 20). 
A rigorous proof of the stability/instability could give 
a linear stability analysis with the corresponding eigen- 
value problem. 

Soliton structures can arise as the result of mod- 
ulational instability of the initial monochromatic UH 
pump wave. We numerically studied the evolution of a 
monochromatic long wavelength UH wave in the frame- 
work of Eqs. (|60|) and (|6ip . Consider the case when a 
monochromatic UH wave of the form 



t=500 



t=1000 



y) = 1.5 exp(0.07ix + 0. 07 iy) 



(68) 



is chosen as an initial condition in a box of side L = 10. 
Thus, the wave Eq. represents the almost homoge- 
neous initial field. The time evolution of the field Eq. 

is shown in Fig. O It is clearly seen the formation 
of one well-shaped soliton, which moves in some direc- 
tion with almost no change in shape and at later times 
coexists with turbulent environment. Depending on the 
amplitude and wave vector of the pump as well as on 
the box length, we observed the formation of two and 
more solitons, and even, in some cases, the formation (si- 
multaneously with solitons) of structures resembling the 
ring-like vortices. No soliton formation, however, was 



(a) >. 10 



(b) 




FIG. 4: (a) Stable evolution of the monopole with A = 1; (b) 
splitting of the dipole with A = 2 into two monopoles. 
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FIG. 5: Evolution of the vortex with m = 1 and A = 0.01 for 
P = 0.25. 



detected in all cases if the amplitude of the initial pump 
was sufficiently small (for example, less than 0.8 for Eq. 
(f68|) ) so that, for a given grid, box size and the wave 
vector of the pump, the threshold of the modulational 
instability was not exceeded. It is important, that in all 
our numerical simulations of the evolution of initial fields 
we did not observed any evidence of collapse. 



VI. CONCLUSION 

In the present paper, we have considered the nonlinear 
interaction between dispersive magnetosonic and high- 
frequency UH waves ion 2D geometry. We have derived 
a set of 2D equations describing the dynamics of non- 
linearly coupled magnetosonic and UH waves. Nonlo- 
cal nonlinearity in the equations prevents collapse and 
results in the possibility of existence of stable 2D non- 
linear structures. We have presented a rigorous proof 
of the absence of collapse in the model under consid- 
eration. We have found numerically different types of 
nonlinear 2D localized structures such as fundamental 
solitons, radially symmetric vortices, nonrotating multi- 
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solitons (dipoles and quadrupoles), and rotating multi- 
solitons (azimuthons). By direct numerical simulations 
we have shown that fundamental solitons are stable and 
do not collapse even with negative initial hamiltonian. 
Multisolitons and vortices decay into fundamental soli- 
tons but can survive over quite considerable times. 



FIG. 6: Evolution of the almost homogeneous initial field 
Eq. (l68|) in the model Eqs. (|60) and (161) . The panels show 
the intensity |$| in the (x, y) plane at times t = 0, t = 1, 
t = 3, and t = 12, respectively. The modulational instability 
of the initial field results in emergence of a self-localized state 
(soliton) . 
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